Bridging the gap: viable reaction pathways from tetrahedrane to benzyne.

The addition of sp-carbon-containing molecules to polycyclic sp3 tetrahedrane (c-C4H4) results in the formation of both o-benzyne (c-C6H4) and benzene (c-C6H6). Since both c-C6H4 and c-C6H6 have been detected in the interstellar medium (ISM), providing additional pathways for their possible astrochemical formation mechanisms can lead to the discovery of other molecules, such as c-C4H4, benzvalyne, and vinylidene (:CCH2). Addition of diatomic carbon (C2), the ethynyl radical (C2H), vinylidene, and acetylene (HCCH) to c-C4H4 is undertaken in individual pathways through high-level quantum chemical computations at the CCSD(T)-F12b/cc-pVTZ-F12 level of theory. The resulting C2 addition pathway proceeds barrierlessly through benzvalyne as an intermediate and reaches a true minimum at c-C6H4, but no leaving groups are produced which is required to dissipate excess energy within an interstellar chemical scheme. Similarly, the C2H addition to c-C4H4 produces benzvalyne as well as its related isomers. This pathway allows for the loss of a hydrogen leaving group to dissipate the resulting energy. Lastly, the HCCH and :CCH2 addition pathways follow through both benzvalene and benzvalyne in order to reach c-C6H6 (benzene) and c-C6H4 (o-benzyne) as well as H2 as the required leaving group. Although there is a barrier to the HCCH addition, the :CCH2 addition presents the contrary with only submerged barriers. These proposed mechanisms provide alternative possibilities for the formation of complex organic molecules in space.

The infrared absorption spectrum of phenylacetylene and its deuterated isotopologue in the mid- to far-IR.

Anharmonicity strongly influences the absorption and emission spectra of polycyclic aromatic hydrocarbon (PAH) molecules. Here, IR-UV ion-dip spectroscopy experiments together with detailed anharmonic computations reveal the presence of fundamental, overtone, as well as 2- and 3-quanta combination band transitions in the far- and mid-infrared absorption spectra of phenylacetylene and its singly deuterated isotopologue. Strong absorption features in the 400-900 cm-1 range originate from CH(D) in-plane and out-of-plane wags and bends, as well as bending motions including the C≡C and CH bonds of the acetylene substituent and the aromatic ring. For phenylacetylene, every absorption feature is assigned either directly or indirectly to a single or multiple vibrational mode(s). The measured spectrum is dense, broad, and structureless in many regions but well characterized by computations. Upon deuteration, large isotopic shifts are observed. At frequencies above 1500 cm-1 for d1-phenylacetylene, a one-to-one match is seen when comparing computations and experiments with all features assigned to combination bands and overtones. The C≡C stretch observed in phenylacetylene is not observed in d1-phenylacetylene due to a computed 40-fold drop in intensity. Overall, a careful treatment of anharmonicity that includes 2- and 3-quanta modes is found to be crucial to understand the rich details of the infrared spectrum of phenylacetylene. Based on these results, it can be expected that such an all-inclusive anharmonic treatment will also be key for unraveling the infrared spectra of PAHs in general.

Interstellar formation of glyceric acid [HOCH2CH(OH)COOH]-The simplest sugar acid.

Glyceric acid [HOCH2CH(OH)COOH]-the simplest sugar acid-represents a key molecule in biochemical processes vital for metabolism in living organisms such as glycolysis. Although critically linked to the origins of life and identified in carbonaceous meteorites with abundances comparable to amino acids, the underlying mechanisms of its formation have remained elusive. Here, we report the very first abiotic synthesis of racemic glyceric acid via the barrierless radical-radical reaction of the hydroxycarbonyl radical (HOCO) with 1,2-dihydroxyethyl (HOCHCH2OH) radical in low-temperature carbon dioxide (CO2) and ethylene glycol (HOCH2CH2OH) ices. Using isomer-selective vacuum ultraviolet photoionization reflectron time-of-flight mass spectrometry, glyceric acid was identified in the gas phase based on the adiabatic ionization energies and isotopic substitution studies. This work reveals the key reaction pathways for glyceric acid synthesis through nonequilibrium reactions from ubiquitous precursor molecules, advancing our fundamental knowledge of the formation pathways of key biorelevant organics-sugar acids-in deep space.

Performance of EOM-CCSD(T)(a)*-Based Quartic Force Fields in Computing Fundamental, Anharmonic Vibrational Frequencies of Molecular Electronically Excited States with Application to the Ã1A″ State of :CCH2 (Vinylidene).

Highly accurate anharmonic vibrational frequencies of electronically excited states are not as easily computed as their ground electronic state counterparts, but recently developed approximate triple excited state methods may be changing that. One emerging excited state method is equation of motion coupled cluster theory at the singles and doubles level with perturbative triples computed via the (a)* formalism, EOMEE-CCSD(T)(a)*. One of the most employed means for the ready computation of vibrational anharmonic frequencies for ground electronic states is second-order vibrational perturbation theory (VPT2), a theory based on quartic force fields (QFFs),fourth-order Taylor series expansions of the potential portion of the internuclear Watson Hamiltonian. The combination of these two is herein benchmarked for its performance for use as a means of computing rovibrational spectra of electronically excited states. Specifically, the EOMEE-CCSD(T)(a)* approach employing a complete basis set extrapolation along with core electron inclusion and relativity (the so-called "CcCR" approach) defining the QFF produces anharmonic fundamental vibrational frequencies within 2.83%, on the average, of reported gas-phase experimentally assigned values for the test set including the A~1A″ states of HCF, HCCl, HSiF, HNO, and HPO. However, some states have exceptional accuracy in the fundamentals, most notably for ν2 of A~1A″ HCCl in which the CcCR QFF value is within 1.8 cm-1 at 927.9 cm-1 (or 0.2%) of the experiment. Additionally, this approach produces rotational constants to, on the absolute average, within 0.41% of available experimental data, showcasing notable accuracy in the computation of rovibronic spectral data. Furthermore, utilizing a hybrid approach composed of harmonic CcCR force constants along with a set of simple EOMEE-CCSD(T)(a)*/aug-cc-pVQZ QFF cubic and quartic force constants is faster than using pure CcCR and better represents those modes that suffer from numerical instability in the anharmonic portion of the QFF, implying that this so-called "CcCR + QZ" QFF approach may be the best for future applications. Finally, complete, rovibrational spectral data are provided for A~1A2 :CCH2, a molecule of potential astrochemical interest, in order to aid in its potential future experimental rovibronic characterization.

Quantum Chemistry and Astrochemistry: A Match Made in the Heavens.

Quantum chemistry can uniquely answer astrochemical questions that no other technique can provide. Computations can be parallelized, automated, and left to run continuously providing exceptional molecular throughput that cannot be done through experimentation. Additionally, the granularity of the individual computations that are required of potential energy surfaces, reaction mechanism pathways, or other quantum chemically derived observables produces a unique mosaic that make up the larger whole. These pieces can be dissected for their individual contributions or evaluated in an ad hoc fashion for each of their roles in generating the larger whole. No other scientific approach is capable of reporting such fine-grained insights. Quantum chemistry also works from a bottom-up approach in providing properties directly from the desired molecule instead of a top-down perspective as required of experiment where molecules have to be linked to observed phenomena. Furthermore, modern quantum chemistry is well within the range of "chemical accuracy" and is approaching "spectroscopic accuracy." As such, the seemingly difficult questions asked by astrochemistry that would not be asked initially for any other application require quantum chemical reference data. While the results of quantum chemical computations are needed to interpret astrochemical observation, modeling, or laboratory experimentation, such hard questions, regardless of the original need to answer them, produce unique solutions. While questions in astrochemistry often require novel developments in and implementations of quantum chemistry as outlined herein, the applications of these solutions will stretch beyond astrochemistry and may yet impact fields much closer to Earth.

A Vision for the Future of Astrochemistry in the Interstellar Medium by 2050.

By 2050, many, but not nearly all, unattributed astronomical spectral features will be conclusively linked to molecular carriers (as opposed to nearly none today in the visible and IR); amino acids will have been observed remotely beyond our solar system; the largest observatories ever constructed on the surface of the Earth or launched beyond it will be operational; high-throughput computation either from brute force or machine learning will provide unprecedented amounts of reference spectral and chemical reaction data; and the chemical fingerprints of the universe delivered by those of us who call ourselves astrochemists will provide astrophysicists with unprecedented resolution for determining how the stars evolve, planets form, and molecules that lead to life originate. Astrochemistry is a relatively young field, but with the entire universe as its playground, the discipline promises to persist as long as telescopic observations are made that require reference data and complementary chemical modeling. While the recent commissionings of the James Webb Space Telescope and Atacama Large Millimeter Array are ushering in the second "golden age" of astrochemistry (with the first being the radio telescopic boom period of the 1970s), this current period of discovery should facilitate unprecedented advances within the next 25 years. Astrochemistry forces the asking of hard questions beyond the physical conditions of our "pale blue dot", and such questions require creative solutions that are influential beyond astrophysics. By 2050, more creative solutions will have been provided, but even more will be needed to answer the continuing question of our astrochemical ignorance.

DFT + F12 QFFs for Cost-Effective Rovibrational Spectral Data Predictions of Ground and Excited Electronic States.

The quest for faster computation of anharmonic vibrational frequencies of both ground and excited electronic states has led to combining coupled cluster theory harmonic force constants with density functional theory cubic and quartic force constants for defining a quartic force field (QFF) utilized in conjunction with vibrational perturbation theory at second order (VPT2). This work shows that explicitly correlated coupled cluster theory at the singles, doubles, and perturbative triples levels [CCSD(T)-F12] provides accurate anharmonic vibrational frequencies and rotational constants when conjoined with any of B3LYP, CAM-B3LYP, BHandHLYP, PBE0, and ωB97XD for roughly one-quarter of the computational time of the CCSD(T)-F12 QFF alone for our test set. As the number of atoms in the molecule increases, however, the anharmonic terms become a greater portion of the QFF, and the cost comparison improves with HOCO+ and formic acid, requiring less than 15 and 10% of the time, respectively. In electronically excited states, PBE0 produces more consistently accurate results. Additionally, as the size of the molecule and, in turn, QFF increase, the cost savings for utilizing such a hybrid approach for both ground- and excited-state computations grows. As such, these methods are promising for predicting accurate rovibrational spectral properties for electronically excited states. In cases where well-behaved potentials for a small selection of targeted excited states are needed, such an approach should reduce the computational cost compared to that of methods requiring semiglobal potential surfaces or variational treatments of the rovibronic Hamiltonian. Such applications include spectral characterization of comets, exoplanets, or any situation in which gas phase molecules are being excited by UV-vis radiation.

Spectroscopy and Photochemistry of OAlNO and Implications for New Metal Chemistry in the Atmosphere.

A new aluminum-bearing species, OAlNO, which has the potential to impact the chemistry of the Earth's upper atmosphere, is characterized via high-level, ab initio, spectroscopic methods. Meteor-ablated aluminum atoms are quickly oxidized to aluminum oxide (AlO) in the mesosphere and lower thermosphere (MLT), where a steady-state layer of AlO then builds up. Concurrent formation of nitric oxide (NO) in the same region of the atmosphere will lead to the bimolecular formation of the OAlNO molecule. Molecular orbital analysis provides fundamental insights into the chemical bonding and energetic arrangement of the triplet (1 3A″) ground state and singlet (1 1A') excited-state species of OAlNO. Additionally, unpaired electrons on the terminal oxygen atom of triplet (1 3A″) OAlNO cause it to be reactive to atmospheric species, potentially impacting climate science and high-altitude chemistry. The triplet (1 3A″) ground-state species exhibits a large permanent dipole moment useful for rotational spectroscopic detection; however, similar rotational constants to the singlet (1 1A') excited-state species will hamper differentiation in a spectrum. Strong infrared intensities will assist in detection and discrimination of the different spin states and isomers. Repulsive electronic excited states of OAlNO will lead to photolysis of the Al-N bond and formation of various electronic states of AlO + NO through nonadiabatic pathways. Reaction through the OAlNO intermediate represents a means for the production of electronically excited AlO, leading to new chemistry in the atmosphere. Excitation to higher-lying electronic states will lead to fluorescence with a minor Stokes shift, useful for laboratory investigation. Such physical properties of this molecule will allow for new, unexplored chemical pathways in the MLT to be considered.

Spectroscopic Constants and Anharmonic Vibrational Frequencies of C(O)OC, c-C2O2 and Their Silicon-Containing Analogues.

Comets are likely to contain various carbon oxide molecules potentially including C(O)OC and c-C2O2 on their surfaces and comae, as well as their silicon-substituted analogues possibly playing a role in the formation of interstellar dust grains. In this work, high-level quantum chemical data are provided to support such potential future astrophysical detection through the generation of predicted rovibrational data. Laboratory-based chemistry would also benefit from such aforementioned computational benchmarking considering these molecules' historic computational and experimental elusiveness. Coupled-cluster singles, doubles, and perturbative triples, the F12b formalism, and the cc-pCVTZ-F12 basis set garner the rapid, yet highly trusted F12-TcCR level of theory leveraged presently. This current work points to all four molecules' strong IR activity, coupled with large intensities, thus suggesting the potential for JWST detection. Although Si(O)OSi possesses a permanent dipole moment significantly larger than those of the other molecules of present interest, the significant abundance of the potential precursor carbon monoxide suggests that the dicarbon dioxide molecules may yet be observable in the microwave region of the electromagnetic spectrum. Thus, this present work details the likely existence and detectability of these four cyclic molecules, providing updated implications compared to previous work performed both experimentally and computationally.

F12+EOM Quartic Force Fields for Rovibrational Predictions of Electronically Excited States.

Quartic force fields (QFFs) constructed using a sum of ground-state CCSD(T)-F12b energies with EOM-CCSD excitation energies are proposed for computation of spectroscopic properties of electronically excited states. This is dubbed the F12+EOM approach and is shown to provide similar accuracy to previous methodologies at lower computational cost. Using explicitly correlated F12 approaches instead of canonical CCSD(T), as in the corresponding (T)+EOM approach, allows for 70-fold improvement in computational time. The mean percent difference between the two methods for anharmonic vibrational frequencies is only 0.10%. A similar approach is also developed herein which accounts for core correlation and scalar relativistic effects, named F12cCR+EOM. The F12+EOM and F12cCR+EOM approaches both match to within 2.5% mean absolute error of experimental fundamental frequencies. These new methods should help in clarifying astronomical spectra by assigning features to vibronic and vibrational transitions of small astromolecules when such data are not available experimentally.

pbqff: Push-Button Quartic Force Fields.

pbqff is an open-source program for fully automating the production of quartic force fields (QFFs) and their corresponding anharmonic spectroscopic data. Rather than being a monolithic piece of code, it consists of several key modules including a generic interface to quantum chemistry codes and, notably, queuing systems; a molecular point group symmetry library; an internal-to-Cartesian coordinate conversion module; a module for the ordinary least-squares fitting of potential energy surfaces; and an improved second-order rotational and vibrational perturbation theory package for asymmetric and symmetric tops that handles type-1 and -2 Fermi resonances, Fermi resonance polyads, and Coriolis resonances. All of these pieces are written in Rust, a modern, safe, and performant programming language that has much to offer for scientific programming. This work introduces pbqff and its surrounding ecosystem, in addition to reporting new anharmonic vibrational data for c-(C)C3H2 and describing how the components of pbqff can be leveraged in other projects.

Electronically excited states of closed-shell azabenzene and purine anion structures: Is 2 D enough for dipole-bound excited states?

Deprotonated azabenzene anions require dipole moments in their corresponding neutral radicals of more than 3.5 D in order to exhibit dipole-bound excited states (DBXSs). This is notably larger than the typical 2.0-2.5 D associated with such behavior. Similar computational analysis on deprotonated purine derivatives also conducted herein only requires the more traditional 2.5 D dipole moment, implying that the single six-membered azabenzene rings have additional factors at play in binding diffuse electrons. The present study also shows that the use of coupled cluster singles and doubles with a double-zeta correlation consistent basis set and additional diffuse functions originating from the center-of-charge for all aspects of the computations decreases the error in predicting DBXSs to less than 0.006 eV at worst and likely less than 0.003 eV for most cases. These results can influence the modeling of molecular spectra beyond fundamental chemical curiosity with application to astrochemistry, solar energy harvesting, and combustion chemistry among others.

Complete, Theoretical Rovibronic Spectral Characterization of the Carbon Monoxide, Water, and Formaldehyde Cations.

New high-level ab initio quartic force field (QFF) methods are explored which provide spectroscopic data for the electronically excited states of the carbon monoxide, water, and formaldehyde cations, sentinel species for expanded, recent cometary spectral analysis. QFFs based on equation-of-motion ionization potential (EOM-IP) with a complete basis set extrapolation and core correlation corrections provide assignment for the fundamental vibrational frequencies of the A˜2B1 and B˜2A1 states of the formaldehyde cation; only three of these frequencies have experimental assignment available. Rotational constants corresponding to these vibrational excitations are also provided for the first time for all electronically excited states of both of these molecules. EOM-IP-CCSDT/CcC computations support tentative re-assignment of the ν1 and ν3 frequencies of the B˜2B2 state of the water cation to approximately 2409.3 cm-1 and 1785.7 cm-1, respectively, due to significant disagreement between experimental assignment and all levels of theory computed herein, as well as work by previous authors. The EOM-IP-CCSDT/CcC QFF achieves agreement to within 12 cm-1 for the fundamental vibrational frequencies of the electronic ground state of the water cation compared to experimental values and to the high-level theoretical benchmarks for variationally-accessible states. Less costly EOM-IP based approaches are also explored using approximate triples coupled cluster methods, as well as electronically excited state QFFs based on EOM-CC3 and the previous (T)+EOM approach. The novel data, including vibrationally corrected rotational constants for all states studied herein, provided by these computations should be useful in clarifying comet evolution or other remote sensing applications in addition to fundamental spectroscopy.

Red-Shifting the Excitation Energy of Carbonic Acid Clusters Via Nonminimum Structures.

Nonminimum carbonic acid clusters provide excitation energies and oscillator strengths in line with observed ice-phase UV absorptions better than traditional optimized minima. This equation-of-motion coupled cluster quantum chemical analysis on carbonic acid monomers and dimers shows that shifts to the dihedral angle for the internal heavy atoms in the monomer produce UV electronic excitations close to 200 nm with oscillator strengths that would produce observable features. This τ(OCOO) dihedral is actually a relatively floppy motion unlike what is often expected for sp2 carbons and can be distorted by 30° away from equilibrium for an energy cost of only 11 kcal/mol. As this dihedral decreases beyond 30°, the excitation energies decrease further. The oscillator strengths do, as well, but only to a point. Hence, the lower-energy distortions of τ(OCOO) are sufficient to produce structures that exhibit excitation energies and oscillator strengths that would red-shift observed spectra of carbonic acid ices away from the highest UV absorption feature at 139 nm. Such data imply that colder temperatures (20 K) in the experimental treatment of carbonic acid ices are freezing these structures out after annealing, whereas the warmer temperature experiments (80 K) are not.

Low-Temperature Thermal Formation of the Cyclic Methylphosphonic Acid Trimer [c-(CH3 PO2 )3 ].

We report the formation of the cyclic methylphosphonic acid trimer [c-(CH3 PO2 )3 ] through condensation reactions during thermal processing of low-temperature methylphosphonic acid samples exploiting photoionization reflectron time-of-flight mass spectrometry (PI-ReTOF-MS) along with electronic structure calculations. Cyclic methylphosphonic acid trimers are formed in the solid state and detected together with its protonated species in the gas phase upon single photon ionization. Our studies provide an understanding of the preparation of phosphorus-bearing potentially prebiotic molecules and the fundamental knowledge of low-temperature phosphorus chemistry in extraterrestrial environments.

Preparation of methanediamine (CH2(NH2)2)-A precursor to nucleobases in the interstellar medium.

Although methanediamine (CH2(NH2)2) has historically been the subject of theoretical scrutiny, it has never been isolated to date. Here, we report the preparation of methanediamine (CH2(NH2)2)-the simplest diamine. Low-temperature interstellar analog ices composed of ammonia and methylamine were exposed to energetic electrons which act as proxies for secondary electrons produced in the track of galactic cosmic rays. These experimental conditions, which simulate the conditions within cold molecular clouds, result in radical formation and initiate aminomethyl (CH2NH2) and amino ([Formula: see text]2) radical chemistry. Exploiting tunable photoionization reflectron time-of-flight mass spectrometry (PI-ReToF-MS) to make isomer-specific assignments, methanediamine was identified in the gas phase upon sublimation, while its isomer methylhydrazine (CH3NHNH2) was not observed. The molecular formula was confirmed to be CH6N2 through the use of isotopically labeled reactants. Methanediamine is the simplest molecule to contain the NCN moiety and could be a vital intermediate in the abiotic formation of heterocyclic and aromatic systems such as nucleobases, which all contain the NCN moiety.

Confirmation of gaseous methanediol from state-of-the-art theoretical rovibrational characterization.

High-level rovibrational characterization of methanediol, the simplest geminal diol, using state-of-the-art, purely ab initio techniques unequivocally confirms previously reported gas phase preparation of this simplest geminal diol in its C2 conformation. The F12-TZ-cCR and F12-DZ-cCR quartic force fields (QFFs) utilized in this work are among the largest coupled cluster-based anharmonic frequencies computed to date, and they match the experimental band origins of the spectral features in the 980-1100 cm-1 range to within 3 cm-1, representing a significant improvement over previous studies. The simulated spectrum also matches the experimental spectrum in the strong Q branch feature and qualitative shape of the 980-1100 cm-1 region. Additionally, the full set of rotational constants, anharmonic vibrational frequencies, and quartic and sextic distortion constants are provided for both the lowest energy C2 conformer as well as the slightly higher Cs conformer. Several vibrational modes have intensities of 60 km mol-1 or higher, facilitating potential astronomical or atmospheric detection of methanediol or further identification in laboratory work especially now that gas phase synthesis of this molecule has been established.

The thermal isomerization of benzvalyne to benzyne.

The isomerization of the highly strained benzvalyne structure to o-benzyne has been investigated using MCSCF and CCSD(T) levels of theory. Two reaction channels were modeled: the disrotatory one which leads directly to the benzyne product, and the conrotatory one which leads to an intermediate which can subsequently lead to the benzyne product. Energies at the MRMP2 level give 22.9 kcal mol-1 for the disrotatory channel and 21.7 and 1.4 kcal mol-1 for the two steps in the conrotatory one. However, the CCSD(T) energies give 19.3 and 14.2 kcal mol-1 for the two conrotatory steps. The first step of the conrotatory channel is significantly higher than the second so is rate determining for this channel. Comparison of the two separate channels shows that the conrotatory has just a slight energetic edge of 1.2 kcal mol-1 at the MRMP2 level. We did not compute the disrotatory channel at the CCSD(T) level due to the significant biradical nature of the wavefunction with natural orbital occupation numbers of 1.2 and 0.8 in the active space.

Theoretical Characterization of Carbonic Acid Clusters in the UV.

Two theoretical structural motifs are proposed to match two experimental solid carbonic acid UV spectra from previous literature ( Astron. Astrophys. 2021, 646, A172): a linear ribbon structure as a single octamer and nonplanar orientations of carbonic acid clusters. The latter have some contribution from approximated amorphous solid carbonic acid in the form of 40 different clusters of 8 carbonic acid molecules ensemble-averaged together, but unoptimized pairs of optimized dimers oriented perpendicular to one another give the strongest intensities of lower energy UV transitions. The linear ribbon structure's predicted spectrum computed with CAM-B3LYP/6-311G(d,p) agrees well with Experimental Solid B─the ß-carbonic acid experimental data in the UV region. Meanwhile, the 40 amorphous clusters are built with a randomization script, and the electronically excited states are calculated with both CAM-B3LYP/6-311G(d,p) and ωB97XD/6-311G(d,p). The resulting theoretical spectrum is constructed by employing a Boltzmann distribution of the intensities and artificially broadening the simulated spectra. The nonplanar dimer pairs are computed with CAM-B3LYP and B3LYP with the 6-311G(d,p) basis set. The results of the amorphous simulation weakly correspond with the Experimental Solid A spectrum, but the fully nonplanar motif matches the experiment much more convincingly. As a result, the previous work appears to have observed the traditional crystalline phase of solid carbonic acid in Experimental Solid B, whereas the nonplanar orientations of the carbonic acids in the clusters appear to correlate with Experimental Solid A. This spectral classification will aid in future laboratory work exploring the role that carbonic acid can play in low temperature, low pressure desorbed environments with potential application to astrochemistry.